ALGEBRAIC SYSTEMS OF MATRICES AND GROBNER BASIS 



GERALD BOURGEOIS August 21,2007 

Abstract. One studies a particular algebraic system where the unknowns 
are matrices. We solve this system according to the parameters values 
thanks to the theory of Grobner basis. 



1. Introduction 

1.1 The problem in question. K is a commutative field of characteristic 
and algebraically closed; n > 2; let J„ and n be the identity and null (n, n) 
matrices. 

Some (n, n) matrices are said ST iff they are simultaneously similar to upper- 
triangular matrices. 

The parameters a, j3, 7 € K are given; one studies the symmetrical system 
in the (n, n) matrices a, b, c with coefficients in K: 

{a + b + c = al n 
a 2 + b 2 + c 2 = 0I n 
a 3 + b 3 + c 3 = jl n 

If a, b, c constitute a solution of S and pairwise commutate, then it's easy to 
prove that r(x) — 6x 3 — 6ax 2 + (3a 2 — 3/3)x + 3af3 — 2j — a 3 is a zero polynomial 
for a,b,c. 

The r's discriminant is (up to the multiplicative factor 216): 

dis{r) = 9a 4 /? - 8a 3 7 - 21a 2 /? 2 + 36a/3 7 - I87 2 - a 6 + 3/3 3 . 
dis(r) = r has a multiple root. 

GL (X 

One carries out in § the changes of functions: a = ai + —I n , b = b\ + —I n , 

O O 

2 

c = c\ + —I n - One deduces that: a\ + b\ + c\ = 0, a\ + b\ + = (/3 ;r)^n 5 

a? + 6? + cf = -J„ where 5 = 2a 3 - 9a/3 + 9 7 . 

Remark 1.1.1. 5 = one of the roots of r is the half sum of the 2 others. 
One can thus rewrite the system as follows: 

{a + b + c = 
a 2 + b 2 + c 2 = al n 
a 3 + b 3 + c 3 = rl n 

with r(x) = 6x 3 — 3<rx — 2r, dis(r) = 3(a 3 — 6r 2 ) and 5 = 9r. 
1.2. The chosen method. In the literature we didn't find any result about 
exact solutions of algebraic system where the unknowns are matrices; on the 
other hand there exist methods of numerical calculation which gives approxi- 
mations of such solutions. 
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Here first one seeks a Grobner basis of the ideal generated by S in the ( non 
commutative) ring of the polynomials in the unknowns a, b, c. For that, one 
uses the formal computation software Bergman (c/. [1]); it will be noted that 
this one works only on homogeneous polynomials. 

Initially one highlights 4 cases according to dis(r) or S is zero or not. Then 
one calculates a Grobner basis in each case. In the continuation, these cases are 
classified from the most simple to the most complex. 

Remark 1.2.1. The non commutative ideals don't admit necessarily finite 

Grobner basis but the ideal generated by § has always a finite basis as it would 

be also the case, for example, for the nonhomogeneous system 

{ab = c,bc = a, ca = b}; on the other hand its associated homogeneous system 

{ab = ct, be = at, ca = bt, at = ta, bt = tb, ct = tc} hasn't any finite Grobner 

basis. 

1.3. Solutions of the system. There exist 4 cases: 

1.3.1. The generic case 6dis(r) ^ 0: the evaluations of r(x) in a, b, c are and 
a, b, c are simultaneously diagonalizable. The resolution of § is thus brought 
back to the case n = 1. 

1.3.2. (5^0 and dis{r) = 0: the evaluations of r(x) in a, b, c are 0; a, b, c pairwise 
commutate and generally aren't diagonalizable; one obtains the general solution 



1.3.3. S = and dis(r) ^ 0: the evaluations of r(x)(x 2 — — ) in a,b,c are 0; 

a, b, c are diagonalizable but generally don't commutate; one obtains the general 
solution of § in 4. 

1.3.4. The " nil-potent" case S = dis(r) = 0: the evaluations of any polynomial 
of degree 5 in a, b, c are null; in general a, 6, c don't commutate but are ST. 

In view of explicit solutions one will use the J.C. Faugere's software FGb 
(c/. [2]) which calculates Grobner basis of ideals in a (commutative) ring of 
polynomials in the unknowns in K. 

There exist solutions such that x 4 isn't a zero polynomial for a, b, c iff n > 9. 

The algebra < a, b, c >, generated by a, 6, c, is a nilpotent semi group; if 
a 4 7^ then its nilpotency class is 5; one exposes some results about the associ- 
ated flag: in particular if n = 9, there exists, up to isomorphism, only one flag 
of length 5. The algebras < a,b,c> are classified up to isomorphism. 
1.4. Special fields. In 6 one proves that if a, j3, 7 are real, if n = 3 and if r 
has only one real root then S hasn't any real solution . 

Then one studies the system 



where v is a known quaternion which isn't real, and a, 6, c are unknown quater- 
nions. One obtains a necessary and sufficient condition on v so that there exist 
solutions which don't commutate with v. To do that one uses software SALSA 
( cf. 4); this software, which is a prolongation of the FGb software, works on 
the real solutions of algebraic systems and its results are certified. 



§ in 3. 
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1.5. Generalization. In 7 one expounds some ideas about the system in the 
matrices (n, n) a, b, c: £ {a k + b k + c k = ctkl n , k = 1, 2, 3, 4}. 

It will be shown that the solutions of the system £ have a more complex 
structure than those of §. 



I s * case: v is invertible and a, b, c commutate with u 2 , v 3 . 

The Bergman software provides a Grobner basis from which one extracts the 
relations: 2.1.0: 2b 2 + 2ab+2a 2 - u 2 = 0, 2.1.1: ab - ba = 0, 
2.1.2: 6a 3 - 3au 2 - 2v 3 = 0. 

Remark 2.1.1. Let v — I n ; then for each a's eigenvalues A, there exists an it's 

eigenvalue n such that 6 A 3 — 3A/x 2 — 2 = 0; it admits a multiple root if /i 6 =6. 

2 nd case: a, b, c commutate only with v 3 . 

This time we obtain relations of which those which follow: 

2.1.3: 2{a 2 + b 2 ) + ab + ba = u 2 , 2.1.4: b 3 a - ab 3 = ba 3 - a 3 b, 

2.1.5: ba 3 - a 3 b = a 2 b 2 + abab - b 2 a 2 - baba. 

2.2 Proposition 1. Let the system in a,b,c 



where u 2 is a known (n, n) matrix which isn't an homothety. 

i) There always exists a solution which doesn't commutate with u 2 . 

ii) If u 2 isn't diagonalizable, then there doesn't exist necessarily some solution 
(a, b, c) of §J such that u commutates with o, b, c. For n ~ 2 a necessary and 
sufficient condition of existence of such a solution is that 6/2 — u 6 isn't nilpotent 
of order 2. 

Proof, i) It's enough to prove it for n = 2; one is reduced to the two following 

cases: 



The FGb software says that §^ always has solutions satisfying u 2 a 7^ au 2 . How- 
ever these solutions don't seem structured. 

ii) I s * case: If u 2 is diagonalizable then one can find a solution such that a, b, c 
are simultaneously diagonalizable with u 2 . 



The matricial polynomial 6X 3 — 'iXu 2 — 2 has as (matricial) resolvant poly- 
nomial q(X) =X 2 -9X + ^-u 6 . 



2. The generic case 5^0 and dis{r) 7^ 
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In the continuation we assume n = 2. 

2 nd case: Let u be such that q isn't a zero polynomial for a matrix which com- 
mutates with u, ie such that the g's matricial discriminant: 6/2 — u 6 isn't a 
square. Then 6/2 — u 6 is nilpotent of order 2 and one can be reduced to 



u = 



e A fi 

„ where e = 6; let a, b, c be a solution which commutates with u; 
el 



3 \ J is the unique 
solution for a. Now one seeks b satisfying 2.1.1 and 2.1.0. Then necessarily 

y 

v o 



'k r 

a is in the form: ( ^ ^ ) and satisfies 2.1.2; 



2 2 / \ 

the discriminant 2u — 3a — \ 3 must be a square what isn't. 



3rd 

case: If u is such that u 2 isn't diagonalizable and I2 w 6 is the square 

6 

of a matrix w, then necessarily spectrum(u) — {a, a} where a ^ 0, u 6 isn't 

diagonalizable, spectrum(l2 — g" 6 ) = {1 — "g - ; 1 — "g - } with a 6 7^ 6; w isn't 

diagonalizable and can be selected as a polynomial in u, spectrum(l2 + w) = 

r a 6 r a 6 
{1 + €\J (1 — ), 1 + eyj (1 — )} where e = ±1; thus I2 + w isn't nilpotent 

and is invertible; it's the cube of an invertible matrix t what one can choose as 
a polynomial in u. Then (Cardan's formula) a = 6~t + 6~u 2 f _1 is a solution 
of 2.1.2 (if w and t are correctly selected) which commutates with u and such 
that spectrum(a) = {A, A} where 6A 3 — 3Aa 2 — 2 = 0. 

Existence of b: if 2u 2 — 3a 2 is nilpotent then 2a 2 = 3A 2 ; the two relations be- 
tween a and A imply a e — 6, which is contradictory. From where we deduce 
the existence of a solution which commutates with u.D 

2.3. Resolution of S. Theorem 1. If S ^ and dis(r) ^ then any solution 
(a,b,c) of S is such that the evaluations of r(x) in a,b,c are and a, b, c are 
simultaneously diagonalizable. 

Proof : according to relations 2.1.1 and 2.1.2, a, b, c pairwise commutate and r 
is a zero polynomial for a, b, c; dis(r) ^ allows to conclude. □ 



3. The case 5 7^ and dis(r) = 

3.1 Installation. Here tr 3 = 6r 2 and t / 0; using a homothety on a,b,c , 

one is reduced to r = 1 and er 3 = 6; after possible multiplication of a, b, c by 

j = e~ or j 2 one may assume a = s/6 and r = 1. 

One placed oneself in C but the reasoning is valid in K. 

Thus all the cases are brought back to a case chosen in advance. 

Here the model is the system obtained for a = = 7 = 1 and that we note §* . 
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3.2. Solution of S*. Theorem 2. If S ^ and dis(r) = then the resolution 
of § is brought back to that of S*; any solution of S* is such that a, b, c are 
simultaneously similar to the block-matrices: 

(I v 0\ h \ /-7 \ 

a' = a J , 6' = 7/, , c' = -a , 
V 0/3/ V° -/V \° V 

where a 2 = 0^,, /3 2 = 0#, 7 2 = V and (p + ip + 8 = n. 

Proof: Relations 2.1.1, 2.1.2 remain valid: thus a,b,c commutate 2 by 2 and 
x 3 — x 2 is a zero polynomial for a, b, c. a 2 ,b 2 7 c 2 are three projectors whose sum 
is 7 n ; by a basis change one is reduced to: 

/% 0\ /0 0\ /0 0\ 

a 2 = J , 6 2 = , c 2 = 1, where tp+ip+6 = n. 

\0 0/ \0 0/ \0 I e J 

Thus ker(a 2 ) = ker(b 2 - I n ) ® /cer(c 2 - I n ) = ker(b - I n ) ker(c - I n ) and 
K n = ker(a - I n ) ® ker(b - 7„) ker(c - I n ) = E @ F G; a, 6, c pairwise 
commutate thus E 1 , F, G are stable by a, b, c; in a basis adapted to this decom- 
position: 

/% 0\ /6i b 2 \ lex c 2 0\ 

a = ai a2 I , 6 = I 7^, , c = c 3 c 4 ; but c = 7„ — a — b 
\0 a 3 04/ V>3 64/ \0 7 ey / 

implies that a, 6, c are as announced in the statement. □ 
4. The case 5 = and dis(r) 7^ 

4.1 Grobner basis. 

!a + b + c = 
a 2 + b 2 + c 2 = u 2 
a 3 + b 3 + c 3 = 

One is interested in the solutions which commutate with u. 

The Bergman software provides a Grobner basis of S2 from which one extracts 

the relations: 4.1.1 a 2 b — ba 2 = 0, 

4.1.2 -bab - a 2 b + 2a 3 - u 2 a = 0, 4.1.3 6a 5 - 5u 2 a 3 + u 4 a = 0. 

4.2 Resolution of the system. Theorem 3. If 5 = and dis(r) 7^ then 

every solution a, b, c of S is such that r(x) (x 2 — — ) is a zero polynomial for a, b, c 

and a, b, c are diagonalizable but generally they don't commutate. 
Moreover 1™ = E F where E, F are a, b, c-stable and such that: 

i) The 7i"s dimension m is even and the restrictions of a, b, c to E are 



simultaneously similar to: A = 




5 






~2~ 2 

Note that these matrices aren't ST. 

ii) r is a zero polynomial for the restrictions of a, 6, c to F and the later are 
simultaneously diagonalizable. 

Proof: Here t = 0, er 7^ and r(x) = 3a;(2x 2 — a). 4.1.3 and cr 7^ imply 
K n = E@F where E = ker(a 2 — —I„) and F = ker(r(a)); one deduces the first 
claim. 

Proof of i): by 4.1.1 E is 6, c-stable; let u £ E be a & 2 's eigenvector: 6 2 (u) = Am 
where A e {-, -,0}. 

c 2 (u) = ou — x 2 (u) — y 2 {u) = — X)u where — A e ^, 0}; this implies 

that A = — and if v e E then 6 2 (v) = c 2 (v) = —v, 

ker(a 2 — —I n ) C ker(b 2 — —I n ) hence 

4.2.1 fcer(a 2 - |j„) - fcer(6 2 - |/„) - fcer(c 2 - |/„). 

By 4.1.2.1 a', 6', c', the restrictions of a, b, c to £7, satisfy: 

-b'a'b' - a' 2 6' + 2a' 3 - era 1 = -b'a'b' - -b' - -a' = 0; thus 

3 3 

= trace(-b' 2 a' + -c') = -trace(-a' + d). 

Then trace(a') — trace{b') — trace(c') — and m is even. 
Thanks to a homothety one may assume that a — 3. 

Thanks to a change of basis we may assume that 0! — diag{Im, we see k 

b' in the following form: 

2 1 T , „ ) ; the conditions b' 2 = (a +b') 2 = I m imply 

V 5 / ¥ + Z) 

3 

Y = Z = m and C/V = --Z^; thanks to a basis change matrix in the form 

4 2 



P 
Q 



( leaving invariant a) we may assume that 



. PUQ- 1 ] ^3 

b' = \ 2 2 1 ; there exist P, Q such that PUQ^ 1 = ^—In 

Qvp- 1 i% " 2 : 



2 



therefore 6' 




These matrices aren't ST because (a'b 1 — b'a') 2 = —3I m . 

Proof of ii): Thanks to a homothety we may assume that a = 2. 
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F = ker(a) ker(a 2 - I n ) = ker(a 2 ) ker(a 2 - I n ). By 4.1.1 F is 6, c-stable. 
Let a" , 6" , c" be the restrictions of a, b, c to F. 

By 4.2.1 r(x) = 6x(x 2 — 1) is a zero polynomial for a",6",c" and the later are 
thus diagonalizable. 

We can conclude thanks to the fact that a" 2 ,6" 2 ,c" 2 are projectors or more 
quickly by the search of a Grobner basis of the ideal generated by the following 
system in a, 6, c where u is a known invertible matrix: 

{a+ b+ c = 
a 2 + b 2 + c 2 = 2u 2 
a 3 - au 2 = b 3 - bu 2 = c 3 - cu 2 = 

We become interested in the solutions which commutate with u; the Bergman 
package provides a basis whose one element is ab — ba = 0. Thus a" , 6" , c" 
commutate 2 by 2 and are simultaneously diagonalizable. □ 

5. The case 6 = and dis(r) = 

5.1 Grobner basis. Let the system 

{a+b+c= 
a 2 + b 2 + c 2 = 
a 3 + b 3 + c 3 = 

The Bergman package provides a sytem equivalent to S4: 5.1.0 a + b + c = 0, 
5.1.1 ab + ba = -2a 2 - 2b 2 , 5.1.2 a 2 b = ba 2 , 5.1.3 2a 3 = a 2 b + bab, 

5.1.4 a 2 ba = a 3 b = -^a 4 , 5.1.5 a 5 = 0. 

5.2 Resolution of §4. Let Jk = [{Jk)ij] be the Jordan nilpotent matrix of 
dimension k {(Jk)ij = except (Jfc)i,i+i = !)• 

Theorem 4. If S = and dis(r) = then every solution a, 6, c of § satisfies: 

5 

i) a 4 = 6 4 ; a6a6 = 6a6a = 2° 4 ' ^ ne va l ue of the other monomials of degree 4 in 
(a, 6) is ~^a 4 - 

ii) a; 5 is a zero polynomial for a, b, c and more generally every monomial of de- 
gree 5 in (a, b, c) is zero. 

iii) There exists a solution such that a, b, c don't commutate 2 by 2 iff n > 3. 

iv) There exists a solution such that a 4 7^ iff n > 9. 

v) < a, 6, c >, the algebra generated by a,b,c, is a nilpotent semi-group and 
a, b, c are ST; if a 4 ^ then its nilpotency class is 5. 

To this semi-group we can link a flag; for n — 9 there exists, up to isomorphism, 
only one such flag of length 5. 

Remark 5.2.1. It's the single case where one doesn't provide explicitly the 
general solution. 

1 5 
Proof: i) 5.1.3 and 5.1.4 imply 2a 4 = a 3 b+abab = —-a 4 +abab and abab = ~^ aA \ 

5 

thanks to the (a <-» b) exchange: baba = -6 4 . 5.1.3 and 5.1.4 imply 
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1 5 
2a 4 = a 2 ba + baba = + baba and baba — -a 4 . 

The remaining relations can be easily proved. 

ii) It's sufficient to prove the result for the monomials in (a, b) because c = —a—b. 
By i) it's sufficient to prove that a 5 = 6 5 = 0; the first is 5.1.5 and the (a «-> b) 
exchange provides the second. 

iii) If n — 2 then the solutions are in the form: a = aJ2, b = PJ2, c = —a — b 
where a, ft € K. 



If n = 3 then a = J3, 6 - 




—a — b ( where x ^ — and 
v r 2 



is a solution whose elements don't commutate, although 
r(x) = x 3 is a zero polynomial for a, b, c. 

iv) The selected method: we choose a Jordan form of a and we seek the matrix 
b; the unknowns are the n 2 b's entries. We rewrite the equations 5.1.1 to 5.1.5 
by interchanging the roles of a and b; we simplify the system with the linear 



relations: 



"ba 



-ba 2 = aba 3 



ba 4 



0; we seek a Grobner basis of 



the ideal generated by this system in the (commutative) polynomial ring in n 2 
unknowns in K. To do that we use the FGb package of J.C. Faugere which is 
consistent with Maple 10. 

a 4 implies n > 5; we review the cases from n = 5 to n = 9. 

n = 5 : a = J5 doesn't provide a solution for b. 

n = 6 : a = diag(J^, 0) doesn't provide a solution for b. 

n = 7 : a = diag(J 5 , J 2 ) and a = diag(J$, O2) don't provide a solution for b. 

n = 8 : a = diag(J 5 , J 3 ), a = diag(J$, J2, 0) and a — diag(J 5 , 03) don't provide 

a solution for b. 

n = 9 : a = diag(J 5 , J4), a = diag(J§, J2, O2), a = diag(Js, O4) and 
a = diag(J$, J 2 , J 2) don't provide a solution for b. 

a = diag(J5, J3, 0) is the only choice which provides solutions; here is a solution 
over Q: 



f° 


1 











3 








°1 




~2 


-1 






4 


















-9 




























1 














X 












-1 








3 





































-1 
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-1 



















































-1 
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-1 














\o 
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The b's value relies on the choice of 9 parameters in K and of 2 parameters 
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in K*. 

v) The first claim is a consequence of ii) ,i) and of the fact that the elements of 
a nilpotent semi-group are ST. 
Now the definitions are those of [3]. 

Definition: a flag T of length I is a filter of I" subspaces: (K)o<i<i 
such that {0} = Vq ^ ... ^ Vi = K™. The JF's signature is the sequence: 
{dim(V i+1 /V))i. 

To a nilpotent semi-group S of nilpotency class I we can associate a flag of 
length I as follows: let {0} c [S^QK")] c ... C [5(K")] C K n where [S*(K n )] 
is the vector space generated by the images of the products Si...Sj where Vj < i 

8j e S. 

If a 4 7^ then the flag associated to < a, 6, c > is of length 5 and by i) 
Vi = a 4 (K n ); moreover Vi < 4 V = [a(V i+1 ,b(V i+1 )]. 

Now we look over the particular case where n = 9 and the flag is of length 5: 

We have seen in iv) that necessarily a ( therefore also 6 and c because 
a 4 = 6 4 = c 4 ) is similar to diag{J§, J3, 0); with Maple 10 we prove that if 
a = diag(J$, J3, 0) then the flag associated to < a, b, c > doesn't depend on the 
6's choice and is: V\ = [ei], V 2 = [ei,e 2 , e 6 ],V 3 = [ei, e 2 ,e 3 , e 6 , e 7 , eg], 
V4 = [ei, e 2 , e 3 , e4, e 6 , 67, e 8 , eg] where {ei, eg} is the K n canonical basis. 
Thus (1,2,3,2,1) is the common signature of the flags of length 5 if n — 9. 

From the precedent calculation we deduce that a and b become upper trian- 
gular in the basis {ei, e 2 , e§, es, ei, eg, e^, es, e^}. □ 

Remark 5.2.2. If a and b are similar, it doesn't imply that c is similar to a as 
this counterexample proves it: let n = 4, a = J4; then (up to order) b is similar 
to a, and c is similar to diag(J 2 , J 2 ). 

Remark 5.2.3. The fact that some solutions o, 6, c are found in the same fixed 
similarity class doesn't imply that linked flags have the same signature as we 
can see if we choose n = 4: to the diag(J2, J 2 ) similarity class we can link flags 
of signature (2, 2) or (1, 2, 1). 

5.3. The algebra < a,b,c >. It's a nilpotent algebra over K whose 
{a, b, ab, ba, a 2 , aba, ab 2 , bab, a 2 b, a 4 } is a vectorial generator; therefore its dimen- 
sion is at most 10. 

Again n = 9 and a = diag(J$, J3,0): then {a, b, ab, ba, a 2 , aba, ab 2 , a 4 } is 
a < a,b,c >'s basis and this algebra is of dimension 8; moreover its cen- 
ter has {a 2 , b 2 , aba, ab 2 , a 4 } as a basis and is of dimension 5. Yet the alge- 
bras obtained according to the choice of b = [bij] aren't isomorphic because 

bab = —aba — Aab 2 + voa 4 and a 2 b = ab 2 — 7^°^ where w = 366,8 — ^2,4; these 
algebras are characterized, up to isomorphism, by the tzj's values. The 6's value 
chosen in the theorem 4 iv)'s proof was obtained for vo = 0. 

6. Special fields 

6.1 The real matrices case. Here we prove that the real variant of S may 
have no solution. 
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Proposition 2. Let a,/3, 7 be known reals; we consider the system: 

{a + b + c — al n 
a 2 + b 2 + c 2 = 0I n 
a 3 + b 3 + c 3 = 74 

where a, 6, c are unknown (n, n) real matrices. 
Assume that the polynomial r has only one real root: 
then §r has at least one solution iff n is even. 

Proof: r has one real root u and two non real roots v ± iw; then dis(r) 7^ 0. 
a) n is even. 



If n = 2 then a solution is a = 11I2, b = 




If n = 2m then a solution is the tensorial product by I m of the former expres- 
sions. 

Remark 6.1.1. More generally if n — 2 and if the 2 nd member of §r is com- 

(v z\ 
z y j ' ^' Z ^ a i somor phic to C, then 

there exists a solution (a, 6, c) € C/ 3 which is unique, up to order, 
b) n is odd. 

After a possible change of the unknown matrices s in — r (s — ul 3 ) 

\(v - u) + iw\ 

we may assume that u = and \v + iw\ = 1. 
I s * case: 8^0 that is v(8v 2 - 9) 7^ or ( because |u| < 1) v ^ 0: 
According to 2. theorem 1, a,b,c are simultaneously diagonalizable over C. 
2 nd case: S = that is 0,±i are the r's roots. This is the case 1.3.3 where 
o = — 2 and r = 0: 

According to 4. theorem 3, F — ker(r{a)) is a C-space vector of odd dimension 
and is a, b, c-stable; moreover the restrictions of a, b, c to F are simultaneously 
diagonalizable over C. 

We can also see F as a M-space vector of odd dimension which is a, b, c stable. 

The following lemma proves that the both cases are impossible: 
Lemma 6.1.2. Let u be a real number, v be a non real complex number and w 
be its conjugate. Let a, b, c be real (n, n) matrices and p be an invertible complex 
matrix such that a = p diag(ai, ■ ■ ■ , a n ) p^ 1 , b = p diag(bi, • • • , b„) p^ 1 , c = 
p diag{c\, ■ ■ ■ ,c n ) p~ x where, for all i, (a,,6j,Cj) is a permutation of (u, v,w). 
Then necessarily n is even. 

Proof: (due to R. Israel) Let U — ker(a — ul n ), U\ = {x £ U;bx = vx}, U2 = 
{x e U;bx = wx}; U = U\ © Ui- The complex conjugation operator leaves U 
invariant and interchanges U\ and Ui. Thus U\ and U2 must have the same 
dimension over R, and therefore also over C. So the dimension of U over C 
is even. Similarly ker(b — ul n ) and ker(c — ul n ) have even dimensions. Since 
C™ = ker(a — ul n ) © ker(b — ul n )®ker{c — ul n ), its dimension n is also even. □ 
6.2 The quaternionic case . Let Ti be the set of quaternions; TC is a real 
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subspace, of basis {l,i,j,k}, of the (2,2) complex matrices. We will write a 
quaternion in the form: x = x\ + x^i + Xsj + x^k. 

Let u, v, w be 3 non all real quaternions which commutate; then there exists 
p G S 2 (p 2 = — 1) such that u, v, w are in the form A + pp where A, /1 e R. 



where the unknown a, 6, c are 3 quaternions. 

If we are interested only by the solutions which commutate with p then, thanks 
to the change of p into i £ C, we return to the system S with K = C and n = 1: 
if the r's roots are (\k+ipk)k<3 then the only solution (up to order) in Ji which 
commutates with p is (A& + ppk)k<3- 

Do there exist solutions in Tt which don't commutate with p ? 

We consider the following example: {u = 0, v E H, w = 1}. 
1 st case: Assume that v <E R; as seen in 2.1. the elements of a matricial solution 
commutate 2 by 2; if (A& + ipk)k<3 are the roots of r(x) = 6x 3 — 3vx — 2, then 
the solutions of Sn are (Afe + ppk)k<3 where p is arbitrary in S 2 . 
If v < \^6, then r has 2 non real roots and E>n has an infinity of solutions. 
2 nd case: We assume in the continuation that: {u = 0, v e H \ M, w = 1}. 

There exists /i G H* such that h~ 1 vh eC\l; thus we may suppose that 
v = vi + iv 2 E C \ R. 

Remark 6.2.1. We may assume that v-i > because if (a, b, c) is a solution 
associated to v, then (a,b,c) is a solution associated to v. 

We have seen that the solutions in H whose the elements commutate with 
v are in C. It remains to see if there exist some solutions in H 3 \ C 3 . We 
calculate a Grobner basis associated to the system whose the 8 real unknowns 
are (aj)j<4, (&i)«<4 with 83 ^ or 04 ^ 0, and the real parameters are (t>i, V2). 
Remark 6.2.2. If we hold a solution, the changes of (03,63) into (—03,-63) 
or of (04,64) into (—04,-64) provide other solutions, ie four solutions for one 
(01, 6i)'s value. 

In the following we use the Salsa package; considering an algebraic system 
over R depending on 2 parameters, it makes it possible to find a partition of 
the parameters plane such that the number of solutions of the system in each 
region is a constant; this number is then calculated by Salsa. It's because the 
number of parameters is bounded by 2 in Salsa that we don't consider a more 
general second member in Sn- 

Three among the basis elements are composed of polynomials in (ai, 61, v\, V2); 
we learn about this subsystem T in the unknowns (ai, 61) with help of Salsa. The 
condition upon v so that T has some solutions is A = 3vi+4i>2 + 18i>ii>| — 18 > 0. 

If A > then there exist six solutions in (ai, 61) which are associated to the 
permutations of a solution in (01, 61, Ci). 

However the fact that T has some solutions doesn't imply the existence 
of S-ft' solutions which don't commutate with v; there exists a second separa- 
tor curve associated to other elements of the basis; the condition to fulfil is 
v\ - 3vf < 0. 



Now we consider the system : 
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For example, if v\ — —4 then has, up to order, 4 solutions which don't 
commutate iff I < v 2 < 4\/3 with I w 2.29. 

What is the number of solutions when v is upon a separator curve ? ( Do 
there exist 2 or more such solutions ?). 



We consider the system where the unknowns are the (n, n) matrices a, 6, c, d 
and where (ak)k<4 C K is known: ^ { flfe + ^ + cfe + d k — a k I n , k = 1,2,3, 4}. 
Now we exhibit 2 cases in which the solutions of S and S are very different. 

7.1. The generic case. Contrary to system S, the calculation of a Grobner 
basis of don't provide any easily manageable relations, and in particular any 
zero polynomial for a. Remember that in the case of S, a, b, c are simultaneously 
diagonalizable ( cf Theorem 1). 

As in 1.1. we can associate to the commutative case a polynomial r, the 
form of which may be turn into r(x) = (x 2 + ux + v)(x 2 — ux + w) where 
u, v, w e K. If x 2 + ux + v ( with roots n, T2) is a zero polynomial for a, c and 
a + c= —ul n and if x 2 — ux + w ( with roots r%, n) is a zero polynomial for 6, d 
and b + d = ul n then a, b, c, d is a solution of ^2, the elements of which aren't 
generally ST. 

Indeed, if n — 2 then we may choose: a — diag(r±, r2), c = diag{r2, r\), 
b = p~ 1 diag(r3,r4)p, d = p~ 1 diag(r4, r-$)p where p is invertible. In the generic 
case we may choose p such that a, b, c, d aren't ST. 

7.2. The "pseudonilpotent" case: a k = 0, k = 1,2,3,4. Remember that, 
in the case of S, a, 6, c are nilpotent matrices ( cf Theorem 4). 
Proposition 3. If n = 2 then the J^'s solutions satisfy one of the two following 
patterns ( up to order): 

i) a, b, c, d are four nilpotent matrices such that a + b + c + d = 0. 

ii) a, b, c, d ( defined up to a multiplicative factor) satisfy: 

7.2.1 {a + b + c + d = 0, a 2 = I 2 , b 2 = jl 2 ( where j = e^r), c 2 = j 2 I 2 , d 2 = 0, 
a + jb + j 2 c = 0}. 

If n > 4 then there exist solutions which don't contain any nilpotent matrix. 
Proof: 1 st case: a is neither nilpotent nor diagonalizable; then we may assume 



2 nd case: a isn't nilpotent but is diagonalizable; then we may assume that 



Remark 7.2.2. For all n if a,b,c,d satisfy 7.2.1 then they constitute a J2' s 
solution where only d is nilpotent. 

With the help of a direct sum we easily deduce the last claim. □ 
7.3. Commentary. The systems § and ^ don't hold similar complexities. 



7. Generalization of S to 4 unknowns 





one of 
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7 bis. A simpler example 



With the help of the Grobner's basis theory, we solve the following system 
where a, b are (n, n) unknown matrices with coefficients in K: 



T 



ab + ba = I n 
ba 2 b = 



Proposition 4. If (a, b) is a solution of T then necessarily n is even (n = 2m) 
and a, 6 are simultaneously similar to 

I m \ (0 q 
z J ' \I m 

where z, g are (m, to) matrices such that zq = qz = 

Proof: ab is a projector then we may assume that ab — diag(I m ,0i) with 

71 

m + l = n; then trace(ab) = to = — ; n is even, Z = to and 6a = diag(O m , I m ). 

I x y\ , In q 
We suppose a = \ , b — ' 



z tj \r s 

where x, y, z, t,p, q, r, s are (to, to) unknown matrices. 

With the Bergman software we study the system in the unknowns x,y,z,t,p,q, 
r, s, u: {xp + yr = u 2 , xq + ys = 0, zp + tr = 0, zq + ts = 0,px + qz = Q,py + qt = 
0, rx + sz = 0, ry + st = u 2 , xu = ux, yu = uy, zu = uz, tu = ut,pu = up, qu = 
uq, ru — ur, su — us}. We obtain (xp + yr)s = 0; here xp + yr = I m and s — 0; 
by the same way x = t = p = 0. We deduce easily that a, b are simultaneously 
similar to 

! ' < l) with zq = qz = and y invertible. 



v z OJ'Xy- 1 
By a change of the /TO(o6)-basis we may assume that y — I m . □ 



8. Conclusion 

Thus we have solved the system § in the 3 first cases and studied, in the 
4 th case, the nilpotent algebra generated by a solution. The use of formal cal- 
culation packages seemed to us necessary in order to obtain such results; in 
particular we mention that the use of Grobner basis in the non commutative 
case is exceptional in the literature ( the Ore algebras case is very different 
because the Ore polynomials satisfy some particular commutation's relations). 
Unfortunately we befall here to the limit of computer's calculation capacity: the 
4 th case with n = 10 imply some overflows of the FGb package; to go further it 
would be necessary to work with a more effective package ( as F5 also devised 
by J.C. Faugere). 

It should be interesting to study in detail the system ^ ( generalization in 
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4 unknowns of S) ; we have seen that the systems § and Y have some solutions 
essentially unalike; moreover Y has Grobner basis too complicated for easy us- 
ing; yes these basis don't triangularize according to the unknowns a, b, c. 

On the other, the problem becomes very complicated if the S's second mem- 
bers aren't homotheties; we have pointed this phenomenon out during the propo- 
sition 1 ( system which has only solutions which don't commutate) and in 6.2. 
where we find requirements so that there exist in the quaternions field some non 
commutating solutions; eventually it should be interesting to study the system 
S-h in the general case. 
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